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The momentum distribution is one of the most important quantities which provides information 
about interactions in many-body systems. At the same time it is a quantity that can easily be 
accessed in experiments on ultracold atoms. In this paper, we consider mixtures of light and heavy 
fermionic atoms in an optical lattice described effectively by the Falicov-Kimball model. Using 
a Monte Carlo method, we study how different ordered density- wave phases can be detected by 
measurement of the momentum distribution of the light atoms. We also demonstrate that ordered 
phases can be seen in Bragg scattering experiments. Our results indicate that the main factor that 
determines the momentum distribution of the light atoms is the trap confinement. On the other 
hand, the pattern formed by the heavy atoms seen in the Bragg scattering experiments is very 
sensitive to the temperature and possibly can be used in low-temperature thermometry. 
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I. INTRODUCTION 

It is well known that the equilibrium momentum distri- 
bution of free noninteracting fermions is described by the 
Fermi-Dirac distribution function. When interactions are 
turned on, the momentum distribution is changed due to 
these interactions. Similarly, if the noninteracting gas is 
confined in a harmonic trap, the distribution can also be 
calculated explicitly [l], however in the limit of a large 
number of fermions it is convenient to use the semiclas- 
sical Thomas-Fermi approximation [2 or other semiclas- 
sical approaches [3 . The situation becomes much more 
complicated if the interaction between fermions cannot 
be neglected. Then, numerical methods usually need to 
be applied, especially when the system is in a trap. On 
the other hand, the momentum distribution function is 
a quantity that can be directly accessed in experiments 
with ultracold atomic gases and therefore it is of a par- 
ticular importance in trying to make contact between 
theory and experiment. 

We examine the problem of a mixture of different mass 
fermionic atoms on an optical lattice at low temperature. 
Since most stable heavy isotopes of alkalies are bosonic, 
we envision mixtures of Li^ with K^^, or light Li^ or 
K^^ mixed with heavy fermionic isotopes of Sr or Yb (it 
turns out that if the heavy particle is a boson with strong 
enough intraspecies repulsion, and one is at a low enough 



temperature, then the statistics of the heavy particle does 
not enter into our model, as it appears effectively like a 
hardcore object). We have already performed numerical 
calculations on this system and shown that one can see 
interesting phenomena similar to viscous fingering as one 
tunes the trap curvature for the light species and moves 
from a phase separated state with the heavies on the out- 
side to one with the heavies on the inside [4 . Interesting 
ordered density- wave patterns appear at low temperature 
in these mixtures. The question is, how do we detect 
the presence of such order with current experimentally 
available techniques? In situations where in situ imag- 
ing with single-site precision is available [5 , one would 
simply look for the ordered phases directly, just as they 
appear in the Monte Carlo snapshots of a particular con- 
figuration of the atoms. But can one see effects of the 
density- wave ordering in a time-of- flight image, or via 
direct Bragg scattering of light off of the density-wave 
pattern? We answer these questions here. 



In Section II, we introduce the model and the tech- 
niques used to solve for the equilibrium properties of the 
mixtures of fermions in a harmonic trap. In Section III, 
we present our numerical results. Conclusions follow in 
Section IV. 
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II. THE MODEL AND METHOD 

A mixture of light and heavy fermionic atoms in a har- 
monic trap (each in one and only one hyperfine atomic 
state and hence acting like a spinless fermion), under the 
assumption that the quantum-mechanical effects of the 
hopping of the heavy atoms can be neglected, is described 
by the Falicov-Kimball Hamiltonian |6l [Tj 

^ = - E (4s- + 4s) + E - 4s 
+ ^(F>-M')«', + f/E4s"^i. (1) 

where c\ is an operator that creates a light fermionic 
atom at site i and U is the on-site interspecies interac- 
tion potential. The symbol Vi is the light particle trap 
potential, and /i is its chemical potential. Vl^ is the trap 
for the heavy atoms and ji^ is its respective chemical po- 
tential. The symbol = or 1 is the number operator 
of the heavy particles, which can be treated as a classical 
variable since the heavy particles do not hop. The same 
Hamiltonian can also describe Fermi-Bose mixtures, pro- 
vided the bosonic atoms are the heavy (localized) ones 
and a strong repulsion between them prevents them from 
occupying one site by two bosons (hard core bosons). 
This model can easily be extended to describe a system 
with soft-core bosons. In such a case Wi = 0, 1, 2, . . . and 
an additional term describing the boson-boson on-site in- 
teraction 



'Hb-B = \U^-B ^ Wi{Wi - 1) 



(2) 



has to be added to the Hamiltonian in Eq. ([T]) [8 , 9 . The 
trap potentials for the light {Vi) and heavy (Vl^) atoms 
are given by 



J 



(3) 



where (xi^yi) is the position of site i {a is the lattice 
constant). The steepness of the potential confining the 
light atoms varies from R = 30a to R = 12.9a, whereas 
it is fixed for the heavy atoms with R^ = 30a. The 
chemical potentials /i and /i^ have been introduced in 
order to control the number of the light and heavy atoms, 
respectively. 

The model is solved by means of a variation of the 
Monte Carlo (MC) method. The method is based on the 
classical Metropolis algorithm modified in such a way 
that systems with both quantum and classical degrees of 
freedom can be simulated [10 . In each MC step, a new 
configuration of the heavy atoms is generated. Then, 
the Hamiltonian in Eq. ([T]) is numerically diagonalized 
to yield all the eigenenergies and eigenstates of the light 
atoms for the trial configuration of the heavy atoms. The 
new configuration of the heavy atoms is accepted accord- 
ing to the same rules as in the Metropolis algorithm, but 



with the free energy of the light atoms used for comparing 
energies instead of the internal energy. The results are 
then averaged over all the configurations generated dur- 
ing the entire MC run. Real-space configurations of the 
atoms for different model parameters (interaction, shape 
of the trapping potential) at different temperatures have 
been studied in Ref. 4, In the present paper we use a 
similar method to investigate the momentum distribu- 
tion and Bragg scattering spectra. Since diagonalization 
of the Hamiltonian in Eq. (TJ for a given configuration 
of the heavy atoms gives all the eigenstates of the light 
atoms, it allows one to also calculate the momentum dis- 
tribution of the light atoms. Under the assumption that 
the optical lattice potential is deep enough that we can 
restrict to a single-band model, then the field operator 
of these atoms ^(r) (expanded in terms of the Bloch 
wavefunctions for the lowest band) is given by 



^(r) = ^Cfc^fc(r), 



(4) 



where = l/N'^.Ciex.p{—ik-Ri) and ^^(r) is the 
Bloch wavefunction. Expanding the Bloch wavefunction 
in terms of the Wannier wavefunctions of the lowest band 
w{r — Ri) which are localized about lattice site Ri yields 
^fc(^) = I/a/^V w{r — Ri) exp {ik • Ri). The summa- 
tion runs over the first Brilliouin zone. Then, the free 
space momentum distribution can be calculated as [TT] 



n{k) = \w{k)\^clcj^, 



(5) 



where w(k) is the Fourier transform of the Wannier func- 
tion. Hence the momentum distribution can be approxi- 
mated by 



n(k) 



N 



k-(Ri-Rj 



(6) 



in the single-band limit, with TV the number of lattice 
sites. 

Here, the quantum mechanical expectation value (...) 
is calculated for a given configuration of the heavy atoms 
and therefore n{k) has to be averaged over the configu- 
rations. 

We ignore dynamic interactions between heavy and 
light atoms due to scattering in the course of the ex- 
pansion in the time-of- flight experiments. This is a rea- 
sonable approximation since the density of the clouds is 
very small once the trap and the optical lattice potentials 
are dropped. Since those potentials are dropped over a 
finite period of time, there can be some effects as the 
potentials are lowered, but these tend to be smaller in 
these systems because the repulsive interspecies interac- 
tion keeps the light and heavy particles away from each 
other in their initial distributions on the lattice prior to 
the time-of- flight experiment. 

Since in our approximation the heavy atoms are lo- 
calized, we are interested in the momentum distribution 
of only the light ones. The heavy atoms, however, may 
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FIG. 1: (Color online) Density profiles of the light (solid red line) and heavy (dashed blue line) atoms for [/ = 5J at a 
temperature k^T = 0.0 IJ. The curvature of the trap for the heavy atoms is given by — 30a, whereas for the light atoms 
it varies from = 12.9a to = 30a. The insets show snapshots of real-space distributions of the heavy (left box) and light 
(right box) atoms. 



display some kind of density-wave ordering, which has 
been demonstrated in Ref. [4] In particular, for some 
regimes of parameters they form a checkerboard pattern 
where the heavy atoms occupy only, e.g.^ black squares. 
In another regime, superpositions of vertical and horizon- 
tal stripes or phase separation, with the heavy and light 
atoms occupying different regions of space have been ob- 
served. Of course, for any nonvanishing interaction be- 
tween both species of atoms the distribution of the light 
atoms is up to some degree correlated with the distribu- 
tion of the heavy ones. As a result, also the light atoms 
are ordered, but unless the interaction is strong enough, 
the magnitude of the density- wave order is much smaller. 
This follows from the fact that the light atoms, in con- 
trast to the heavy ones, gain kinetic energy while delo- 
calized. 

The ordering of the heavy atoms can be analyzed by 
use of scattering of light (Bragg scattering, see e.g.^ Ref. 
[12). This kind of experiment is an equivalent to neutron 
or X-ray diffraction for the solid state. However, due 
to the difference of lattice constants between solid state 
crystals and optical lattices, the required wavelength cor- 
responds to that of visible light. The observation of well- 
defined Bragg peaks has been used to confirm a crys- 
talline structure formed by atoms in an optical lattice 
[12]. Since this method allows one to determine the dis- 
tance between atoms forming the crystalline structure, it 
can be applied in order to detect the checkerboard pat- 
tern as well. In the case of strong repulsive interaction 
between the light and heavy atoms, checkerboard squares 
of different color are occupied by different kinds of atoms 
and the problem is similar to the detection of antiferro- 
magnetic (AF) order. It has recently been proposed to 
use Bragg diffraction of light to detect such an order [13] 
in a Hubbard system. In the case of AF order, the spin- 
dependent scattering is achieved by using the probe light 
frequency near atomic resonance, where the interaction 
between light and atoms is neither purely diffractive nor 
purely absorptive. In a similar way the probe light can be 



tuned to be scattered in a different way by the light and 
heavy atoms. As a result, we can expect a (tt, tt) peak 
in the case of a two-dimensional checkerboard pattern 
with the light scattered by the heavy atoms. Moreover, 
the intensity of this peak can give information about the 
fraction of the system occupied by ordered atoms. This 
method can also be used to detect other types of correla- 
tions. If the "labyrinthine" patterns obtained in Monte 
Carlo simulations [4] are superpositions or mixtures of 
horizontal and vertical stripes, the Bragg scattering in 
this case should reveal (0, tt) and (tt, 0) peaks. 

The Bragg spectra integrated over the frequency gives 
the static structure factor S{k)^ a key quantity that can 
be expressed as 

S{k) = ^^w,e^''-^^. (7) 

i 

In a practical realization we get a new configuration {wi} 
in each MC step and in each step S{k) is calculated. 
Then, S{k) is averaged over the entire MC run. 

III. RESULTS 

Simulations have been carried out on a 50 x 50 square 
lattice with 625 heavy and 625 light atoms. In contrast 
to solid state simulations, there is no translational sym- 
metry in the real system and therefore we use hard-wall 
boundary conditions at the edges. For steep harmonic 
trap potentials (and at low T) most of the atoms are in 
the center of the system, however, if the potentials are 
shallow the results may be affected by the finite size of 
the cluster. Fig. [T] illustrates the real-space distributions 
of the light and heavy atoms for different shapes of the 
trap for the light atoms [see Eq. (|3| for the trap potential 
parametrization] . The trapping potential for the heavy 
atoms is kept at = 30a for all simulations. 

From Fig. [ijone can see that if the light atoms' trap 
potential is sufficiently steep and the repulsion between 
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FIG. 2: (Color online) Momentum distribution of the light atoms for the cases shown in Fig. T] The upper row shows a 
comparison of cross sections along the x axis for U = 5J (solid red line) and noninteracting case (dashed black line). The inset 
in the first plot shows the momentum distribution of the light atoms confined in a circular well with infinite walls (see text). 



the two species of atoms is relatively strong {U = 5J), 
the heavy atoms are pushed out from the center of the 
trap. This is an intuitive result. However, when the 
trap potential for the light atoms becomes shallower they 
start to spread and occupy the periphery of the cluster 
even before the trap curvatures are set to be equal. For 
identical potentials, the heavy atoms are concentrated 
in the central part of the trap, surrounded by the light 
atoms. In such a configuration the light atoms can gain 
kinetic energy at the expense of the potential energy in 
the harmonic trap. 



A. Momentum distribution of the light atoms 

Fig. |2] shows the momentum distribution of the light 
atoms for the same parameters as in Fig. [l] Since the 
strong repulsion leads to phase separation, at least in 
the limiting cases, the light atoms are able to move freely 
within the region not occupied by the heavy atoms. Nev- 
ertheless, one can notice that there is no sign of the Fermi 
surface in the momentum distribution, which results from 
the inhomogeneity of the system. The upper row of pan- 
els presents a comparison of the crossections of the actual 
distribution with that of a noninteracting fermionic gas 
in a harmonic trap. The difference results from the inter- 
action between the light and heavy atoms. Generally, for 
all the shapes of the trap potential, the maximal value of 
the momentum distribution function at k = of the light 



atoms is reduced with respect to the noninteracting case. 
For a relatively steep trap potential for the light atoms, 
the interaction leads to an occurrence of additional fea- 
tures in the momentum distribution [see the red (solid) 
curve in the upper left panel of Fig. |2]. The occurrence 
of these additional inflection points results from a further 
confinement of the light atoms in an area surrounded by 
the heavy atoms, as can be seen in Fig. |3] In order to 




FIG. 3: (Color online) Comparison of the real-space distri- 
bution of light atoms with (left graph, U = 5 J) and without 
(right graph) interaction with heavy atoms (R = 12.9a). 

confirm the origin of these features, we calculated the 
momentum distribution of the light atoms in an infinite 
round well with a diameter equal to the inner diameter 
of the ring formed by the heavy atoms. In the resulting 
distribution, presented in the inset in Fig. [2j the addi- 
tional peak in the center of the distribution is even more 
pronounced. 

Independent of the strength of the interaction, the in- 



homogeneity destroys any signs of the Fermi edge in the 
momentum distribution. When the trap becomes shal- 
lower, behavior that looks like a Fermi edge appears, but 
this is due to the box boundary conditions at the edge of 
the lattice. 

Even though there is no explicit signal in the momen- 
tum distribution function which shows the presence of 
different kinds of density wave ordering, the momentum 
distribution function does have a strong dependence on 
the appearance of phase separation, as can be seen in the 
different distribution functions in Fig|2] For example, it 
is known that in a homogeneous lattice, the momentum 
distribution function must decrease below the noninter- 
acting value for small momentum and increase above the 
noninteracting value for large momentum [14 when the 
system phase separates at low temperature. This be- 
havior is clearly seen in all of the data, and most likely 
is arising from different forms of quantum confinement 
effects associated with the phase separation. Unfortu- 
nately, it is not easy to disentangle this phase separation 
effect from the effect of the trap, which has a similar 
effect on the momentum distribution function, except in 
the case where the heavy particles surround the light ones 
and confine them with a sharp boundary. In that case, 
the phase separation effect causes a "dimple" in the mo- 
mentum distribution function near k = 0. 

If the interspecies repulsion is too weak to lead to phase 
separation, a checkerboard configuration may be formed. 
Fig. [4] shows examples of real-space configurations of the 
light atoms for U = J. It turns out that in this case 
the momentum distribution is hardly affected by the in- 
teraction. Fig. [5] presents a comparison of the momen- 
tum distribution at temperatures above and below the 
temperature, at which a checkerboard pattern is formed. 
When the temperature is lowered, the momentum dis- 
tribution is enhanced for small /c, though the difference 
much smaller. Moreover, formation of regions with den- 
sity wave order does not affect the distribution in a sig- 
nificant way. 



B. Heavy atom configurations 

Since the heavy atoms are localized in the proposed 
approach, one cannot analyze their momentum distribu- 
tion. Instead, we determine structure factors defined by 
Eq. [7[ In particular, we are interested in how different 
density-wave-ordered patterns are reflected in the struc- 
ture factor. It is well known from solid state physics that 
different orderings produce characteristic Bragg spectra. 
However, in the case of cold atoms the spectra are addi- 
tionally affected by the confining potential that keeps the 
atoms inside the trap. Neglecting the interspecies inter- 
action that may lead to phase separation or ordering, at 
low temperature, the heavy atoms occupy the bottom of 
the trap forming a circular region. Fourier transforma- 
tion of such a configuration consists of a finite- width peak 
at position k = (0, 0). Plotted in a region (0, 27r) x (0, 27r), 
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FIG. 4: (Color online) Real-space density distribution of the 
light atoms for [/ = J and /cbT = O.OOOIJ. The upper panel 
(a) corresponds to the trapping potential with i? = 12.9a and 
the lower (b) to R = 17a. 



the peak splits into four quarters which are visible in 
each corner [at k = (0,0), (0,27r), (27r, 0) and (27r,27r)]. 
The peak itself simply results from a finite number of 
heavy atoms, since according to Eq. ([7|), S{0) is equal 
to the average concentration of the heavy atoms. The 
gathering of atoms in the center of the trap broadens 
it, as can be seen in Fig. [6j When temperature or in- 
teraction moves the atoms to more peripheral areas, the 
peak width shrinks, finally taking on a delta-function- 
like form for a random distribution of atoms. It changes 
the spectral weight around k = (0,0) (and equivalent 
points), affecting the magnitude of the peaks represent- 
ing ordered phases. This will be discussed in more detail 
below. Since we are mainly interested in patterns formed 
by the heavy atoms, we will neglect the k = (0, 0) peak 
and focus on the remainder of the spectrum. Neverthe- 
less, even in the presence of ordering, this peak is still 
the most pronounced feature of the spectrum. 



Most of the results have been obtained at finite tem- 
peratures by means of the Monte Carlo method. These 
results can also be compared with zero-temperature lo- 
cal density approximation (LDA) results, which we now 
describe. 
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FIG. 5: (Color online) Momentum distribution of the light atoms for the cases shown in Fig. [4] The upper row shows a 
comparison of cross sections along the x axis for /cbT = O.OOOIJ (dashed black line) and /cbT = 0.005 J (solid red line). The 
insets blows up the momentum distribution for atoms with energies close to the Fermi energy. The momentum distributions 
shown in the lower row were determined at /cbT = lO"'^ J. 




FIG. 6: (Color online) Left panel: ground state real-space 
configuration of the heavy atoms in a harmonic trap. Right 
panel: corresponding structure factor. 



1. Zero temperature results: LDA 

We construct the LDA at T = from the homoge- 
neous, grand canonical phase diagram, where the ground 
state phases are given as a function of the chemical po- 
tentials of heavy and light atoms [15]. The procedure 



is as follows. For each lattice site, we determine the lo- 
cal chemical potentials by subtracting the trap poten- 
tial at that lattice site from a trial global chemical po- 
tential. The global chemical potential is then adjusted 
to produce the correct total number of heavy and light 
atoms in the trap. Next, using the local chemical poten- 
tials, we map out the homogeneous phase diagram for 
each site within the trap. Two sets of model parameters 
that give nontrivial configurations have been analyzed: 
R = 12.9a, U = J said R = 18.5a, U 5J. In the 
former case, a checkerboard-type configuration is formed 
in the center of the trap, where the heavy atoms occupy, 
let us say, the black squares, and the light atoms are pri- 
marily on the white ones. The central part is surrounded 
by rings of different phases. In the latter case, the center 
of the trap is occupied by light atoms, while the heavy 
ones form a relatively thick ring composed mainly of var- 
ious stripe phases. The spatial distributions have been 
obtained taking into account all possible periodic phases 
with unit cells consisting of no more than 4 lattice sites 
(in all possible shapes). The candidate phases are pre- 
sented in Fig. [7| The zero-temperature LDA configu- 
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FIG. 7: (Color online) All the configurations of the heavy atoms that have been taken into account in the LDA calculations. 
The letters correspond to different regions presented in Figs. [8j\ and[ToJ^. 



rations are presented in the lower panels in Figs. [8]A 
and P^OK; where different letters correspond to different 
phases marked with the same identifying letters as in 
Fig. [7j The Bragg spectra presented in the upper panels 
m Figs. ^ and[lO}V are calculated in the following way: 
first, the positions of the 5- function- type peaks are deter- 
mined from the Fourier transforms of the configurations 
depicted in Fig. [7| then the relative spectral weight pro- 
portional to the fraction of the trap occupied by the cor- 
responding phase is assigned to the given peaks. Finally, 
the peaks are slightly broadened to make the presenta- 
tion more clear. This is necessary because the LDA does 
not know about the finite size of the system and hence 
always displays perfect delta-function peaks. 

It can be seen in Fig. [8]A, that for R = 12.9a and 
U = J di checkerboard phase occupies the central part 
of the trap. It leads to a highly pronounced peak at 
k = (7r,7r). Other phases give smaller peaks at (7r,7r/2), 
(7r,0), (27r/3, 27r/3) (and points obtained by symmetry 
operations). 

For R = 18.5a and U = 5J there are no heavy atoms 
in the center of the trap: they are distributed in a ring 
of some width with two different phases that have axial 
stripes on its innermost part. Since the stripes are invari- 
ant with respect to translations along the axis (neglect- 
ing the finite size of the system), all peaks are located at 
kx = OT ky = 0, namely at = (0, tt), (tt, 0), (0, ±27r/3) 
and (±27r/3,0). The outer part of the ring is densely 
filled with heavy atoms and its Fourier transform con- 
tributes to the peak at /c = (0, 0). 



2. Finite temperatures: Monte Carlo 

In order to investigate how temperature affects the pat- 
terns formed by the heavy atoms, Monte Carlo simula- 
tions have been carried out. For a sufficiently long MC 



run, a number of independent heavy atom configurations 
are generated. The length of a single MC run depends on 
the temperature and the model parameters (interaction, 
shape of the harmonic trap), but usually it is on the order 
of 10^ MC steps from which about 10^ — 10^ independent 
configurations have been selected. For each configuration 
the Fourier transform has been calculated. Figs. [8pCD 
and [TO^CD present results averaged over all configura- 
tions generated for a given set of model parameters. We 
need to comment about the averaging procedure. In an 
experiment, the Bragg spectra of a single configuration 
can be observed. However, in many cases the configu- 
ration (and the corresponding spectra) changes signifi- 
cantly between successive "snapshots" . In order to make 
our results independent of any particular distribution of 
the atoms, we decided to present results that can char- 
acterize the system at a given temperature and model 
parameters. Fourier transforms of very similar configu- 
rations often have similar shape, but may have different 
sign, depending on the details, e. g.^ on the phases of the 
order parameter for a checkerboard phase. Therefore, in 
Figs. |8^CD and[TO^CD, we present averaged absolute 
values of the spectra, given by 



S{k) 



1 



MN 



M 

E 



N 



(8) 



where N and M are the number lattice sites and the 
number of MC "snapshots", respectively and Wi^^n is 
equal to or 1, depending on whether site i is occu- 
pied or empty in the m-th "snapshot". Additionally, 
the resulting Bragg spectra are self averaged using ro- 
tational and reflection symmetries of the lattice and the 
trap. This means that the presented spectra are calcu- 
lated as averages of S{kx^ky)^ S{kx^—ky)^ S{—kx^ky), 
S{ Siky^kx)-) S{ky^ kx)-, ky^kx)-) and 

S{—ky^—kx). Since we are not interested in the k = (0, 0) 
peak (and equivalent peaks in the remaining corners). 
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FIG. 8: (Color online) Upper row: Fourier transforms of the spatial distributions of the heavy atoms; Lower row: corresponding 
spatial distributions for U — J. Figure A shows a transformation of a configuration obtained within the LDA. In the figure 
presenting the corresponding real-space configuration regions occupied by different phases are filled with different colors and 
marked by letters. Each letter stands for one configuration and the configurations are shown in Fig. |7| Figures B, C, and D 
correspond to results at finite temperature, /cbT = 0.005J, O.OIJ, and 0.05J, respectively. 
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FIG. 9: (Color online) Temperature dependence of the Bragg 
spectra in the center of the trap. More precisely, this figure 
shows "^'j^ S{k) — S{0, 7v) where the sum is calculated over the 
square marked by white dotted line in Fig. [S] 



the false-color scales in Figs. [8] and [To] are chosen in such 
a way that (independent of the temperature) only the 
peaks resulting from ordering are correctly represented. 
But of course, the scale is kept the same in panels B, C, 
and D. 

In the case presented in Fig. [sjA (LDA, U = J, R = 
12.9a), most of the heavy atoms form a checkerboard 
configuration leading to a strong peak at A; = (7r,7r). In 



the peripheral areas, the concentration of heavy atoms 
increases and they form more dense patterns, which in 
turn, lead to less pronounced peaks. One can notice that 
the central peak in Fig. is fourfold split. This results 
from an imperfection of the checkerboard ordering of the 
heavy atoms. This can be explained from an analysis 
of the LDA: the effective chemical potential varies con- 
tinuously when the distance from the center of the trap 
increases. Hence one expects that the density of atoms 
(both light and heavy) should also vary in a continuous 
way. However, in some regions, the checkerboard order- 
ing minimizes the energy, hence in that ordered region 
the density (at least of the heavy atoms) is constant. A 
competition between these two tendencies leads to flaws 
in the center of the trap. Due to topological reasons, such 
a crack has to run across the entire ordered area. Then 
the question is why the fourfold split is not visible in the 
LDA results? This is connected with the maximum size 
of the unit cell taken into account in the LDA calcula- 
tions. The 4 site unit cell is too small to describe a large 
checkerboard area with a single line defect, but by using 
a larger unit cell, one should see such a configuration. 

As can be expected, the ordering is reduced when the 
temperature increases and so are the corresponding fea- 
tures in the Bragg spectra. Figs. |8p and [8p) illustrate 
how the multi-peak structure is smoothed out by destroy- 
ing the real-space ordering. In order to describe this pro- 
cess in a more quantitative way, we calculated how the 
spectral weight of the central peak(s) decreases with in- 
creasing temperature. Figure |9] shows the spectral weight 
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FIG. 10: (Color online) The same as in Fig. [s] but for U = 5J. Panels B, C, and D correspond to temperature /cbT = 
O.OIJ, 0.05J, and 0.15J, respectively. 
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FIG. 11: (Color online) Temperature dependence of the Bragg 
spectra at the edge of the trap: this figure shows Ylk ^i^) ~ 
S{7t,7t) where the sum is calculated over the rectangle marked 
by red dotted line in Fig. [To] 



in a region (O.TStt < < 1.257r) x (0.75 < ky < 1.257r) 
(marked by the white square in Fig. |8]A.) as a function 
of temperature. The plot is shifted in such a way that 
the weight ai k = (tt, 0) is equal to zero. The shift is 
necessary since with increasing temperature less and less 
weight is attributed to the peak at /c = (7r,0), which in- 
creases the reference level more than should result from 
the reduction of the central peak (this effect is especially 
visible in Fig. 



10, where the color around (7r,7r) becomes 



very bright at high temperature) 



IV. SUMMARY 

In this paper, we have examined two simple experimen- 
tal probes that can reveal information about ordered den- 
sity wave phases in mixtures of light and heavy fermionic 
atoms (or equivalently light fermionic and heavy "hard 
core" bosonic atoms) on an optical lattice. Namely, we 
examined the momentum distribution function, which 
comes from a time-of-fiight expansion experiment and 
the Bragg scattering signal that would come from scat- 
tered optical light that scatters off of the density wave 
pattern. 

The momentum distribution function does not pro- 
vide significant information about various ordered den- 
sity wave phases, but it does provide some intuition about 
phase separated states. In particular, as the distribution 
flattens and broadens, one has an indication of phase 
separation setting in. Furthermore, if the heavy atoms 
surround the lights and confine the light atoms with es- 
sentially a hard wall boundary condition, then the mo- 
mentum distribution function develops a sharp dimple at 
low momentum which does provide a characteristic shape 
signalling that form of phase separation. 

Bragg scattering is much more effective at showing the 
presence of ordered density wave phases, as new Bragg 
reflection "spots" appear at appropriate ordering vectors 
for the different types of order present in the sample. 
The weight underneath these peaks is proportional to 
the strength of the ordering, and to the volume of re- 
gions which are ordered, and hence they can be used for 
accurate low-temperature thermometry of these systems. 
As T is lowered, the weight in the peaks grows and can 
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be calibrated via numerical calculations to produce an 
appropriate temperature of the system. 

One caveat of this work, however, is that one must 
cool the mixture down to a low enough temperature that 
the ordering appears in the system. Typically, this re- 
quires a temperature at least as low as about l/40th of 
the bandwidth, and often substantially lower. This is an 
aggressively low temperature with current cooling tech- 
nology, but hopefully can be reached as it becomes eas- 
ier to manipulate entropy distributions within trapped 
atomic systems. Finally, we also should note that direct 
in situ imaging via apparatus like the quantum gas mi- 
croscope would provide even more convincing pictures of 
the ordered phases, and the fluctuations about that or- 
der, than the above proposed methods, but we are not 



aware of any plans to examine these kinds of mixtures 
with such machines in the near term. 
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